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THE  STRESS  DEPENDENCE  OF  ADSORPTION 

MIKAEL  CIFTAN  and  EDWARD  SAIBEL 

Duke  University,  Durham,  N.  C.  27706  and 

U.  S.  Army  Research  Office,  Research  Triangle  Park,  N.  C.  27709 


ABSTRACT 

Considering  a  gas-solid  adsorption  system,  we  present  a  simple  analysis  of  the 
changes  in  the  chemical  potential  of  the  adsorbed  gas  as  the  substrate  solid  is 
subjected  to  external  forces.  Using  the  canonical  partition  function,  we  show 
that  the  chemical  potential  varies  with  applied  stress  which  variation  we  call 
the  chemostress  coefficient.  Numerical  values  obtained  for  these  coefficients 
indicate  that  the  magnitude  of  the  effect  depends  sensitively  on  the  solid-gas 
interaction  potential  which  in  turn  depends  on  the  changes  of  the  polarizabilities 
and  of  excitation  energies  of  the  solid. 


INTRODUCTION 

It  has  always  been  an  open  question  in  all  branches  of  science  whether  in 
developing  a  theory  to  describe  a  given  phenomenon  one  has  taken  into  account  the 
minimum  number  of  fundamentally  Independent  variables.  Ideally  the  more  number 
of  such  variables  a  theory  is  based  on,  the  more  detailed  description  must  ensue 
from  it  so  that  these  variables  are  not  considered  as  "just  some  parameters  to 
fit  the  data." 

A  remarkable  example  of  discovery  of  Buch  fundamental  variables  is  the  spin  of 
an  electron, and  more  such  fundamental  internal  degrees  of  freedom  have  been  dis¬ 
covered  in  elementary  particle  physics.  On  the  classical  statistical  thermody- 
/  namlc  level  one  has  derived  guidance  from,  for  example,  the  Gibbs  phase  rule. 
However,  when  we  start  considering  non-equilibrium  phenomena  on  a  statistical 
physical  level,  we  find  that  the  classical  phenomenological  equilibrium  thermo¬ 
dynamic  concepts  do  not  suffice  to  give  the  details  that  we  want  to  obtain  from 
the  system  under  study.  In  this  paper  we  shall  take  up  the  concept  of  the  chemi¬ 
cal  potential  that  was  considered  by  Gibbs  but  we  shall  generalise  it  so  that  it 
encompasses  Interactions  that  Involve  coupling  between  mechanical  and  chemical 
energy  contents  of  the  system.  We  shall  see  that  the  ordinary  chemical  potential 
so  generalized  will  depend  on  mechanical  stresses  within  a  solid  or  at  the  gas- 

«*  or  2«  i#a 
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solid  interface. 

The  need  for  such  a  generalization  arises  in  connection  with  a  number  of  pheno¬ 
mena  in  which  in  fact  a  coupling  of  mechanical  and  chemical  effects  have  been  ob-  * 
served.  Among  these  effects  are:  the  Rebinder  effect  which  shows  a  change  in 
microhardness  of  a  solid  when  the  surface  of  the  substance  is  in  contact  with 
certain  classes  of  chemicals;  the  well  known  phenomenon  of  stress-corrosion  which 
shows  highly  increased  rate  of  corrosion  of  a  substance  in  certain  chemical  en¬ 
vironments  and  under  stresses  such  that  either  the  chemical  or  the  stresses  alone 
do  not  produce  a  significant  rate  of  corrosion. 

In  a  series  of  papers  (ref.  1  to  6)  we  have  derived  expressions  for  the  depen¬ 
dence  of  the  chemical  potential  on  stresses  in  solids  for  the  adsorption  phenomenon. 
This  interaction,  which  we  call  the  chemostress  effect,  now  gives  a  level  of  detail 
of  the  dependence  of  this  thermodynamic  quantity  on  internal  variables  of  the  sub¬ 
stances  such  as  polarizabilities  of  the  atom  and  the  substance,  phonon  spectrum 
of  the  solid  substrate,  excitation  spectra,  band  structure.  Interatomic  potential 
characteristics,  variation  of  these  with  pressure,  etc.  Such  details  which,  tie 
the  adsorption  phenomena  to  chemical  and  physical  specificities  of  £he  speeies 
involved  were  not  available  in  previous  theories. 

Furthermore,  a  number  of  specific  new  experiments  that  need  be  carried  out  for 
new  quantities  that  come  out  of  this  theory  are  shown  in  these  papers.  Among 
these  we  mention  a  few  here  such  as  the  change  of  heat  of  adsorption  with  stress, 
extension  of  the  "Jarman  Rule"  of  polarization  from  alkali  halides  to  other  solids, 
variation  of  phonon  and  excitation  Bpectra  with  stress. 

Another  aspect  of  the  theory  is  that  it  becomes  a  fertile  ground  for  testing 
several  other  theories  such  as  the  theory  of  lonlclty,  oxidation  theories,  hydrogen 
embrittlement. 

Because  we  have  discussed  in  previous  publications  (ref.  1-6)  details  of  several 
methods  of  calculating  the  chemostress  coefficient  and  because  need  exists  to 
extract  from  these  publications  the  minimal  theoretical  cops.  tat  still 

supports  the  existence  of  this  chemostress  effect,  we  shall  r  •«..  <  orial  pre- 

t 

sentatlon  of  the  simplest  of  these  theoretical  treatments  bjlow. 

A  comprehensive  review  of  the  statistical  theromodyn-mics  of  adsorption  has 
been  given  by  Steele  (ref.  7).  We  make  extensive  use  of  the  works  of  Steele  and  ^ 
Halsey  (ref.  8),  Freeman  and  Halsey  (ref.  9),  Sams,  Constabaris  and  Halsey  (ref.  10) 
and  of  some  of  the  more  general  findings  of  Steele  (ref.  11).  The  present  study 
refers  not  to  the  dependence  of  the  chealcal  potential  v  on  the  pressure  of  the 
adsorbed  gas  but  on  the  stresses  sppllad  to  the  solid  substrate.  Nor  does  the 
present  study  concentrate  on  the  specificities  that  would  be  associated  with 
particular  gas-substrate  systems,  specificities  that  would  result  from  details 
of  the  interatomic  potential  peculiar  for  each  system,  the  crystallographic 
orientation  of  the  surface,  degree  of  Imperfections  within  the  solid  and  on  the 
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surface,  etc.  Such  a  detailed  study  of  special  classes  of  metals,  ionic  solids, 
etc.  with  specified  imperfections  is  deferred  to  the  future.  Rather,  the  signi¬ 
ficant  point  of  the  present  analysis  is  the  result  that  even  with  a  rather  idealized 
solid  or  a  non-metal,  the  chemical  potential  of  the  adsorbed  gas  may  vary  1-101 
with  moderately  high  pressures  of  the  order  of  3-8  kilobars  applied  to  the  solid. 

It  remains  to  be  seen  whether  imperfections  in  the  solid  increase  or  decrease  the 

magnitude  of  this  chemostress  effect.  We  argue  that  within  the  limits  of  this 
0  * 

uncertainty,  this  coefficient,  or  rather  this  tensor  (since  it  depends  on  crystal¬ 
line  surface  orientation  and  the  tensorlal  nature  of  the  stresses)  may  be  a  non- 
negligible  quantity.  We  also  emphasize  that  the  chemostress  effect  need  not  be 
confined  to  physisorption  and  chemisorption  processes  alone,  that  it  is  a  more 
general  effect. 


THE  CHEMICAL  POTENTIAL  FOR  ADSORBED  GAS 

Starting  with  the  canonical  ensemble,  expressions  for  the  variation  of  the 
chemical  potential  with  stress,  which  we  shall  call  the  chemostress  coefficient, 
will  be  derived.  The  numerical  results  obtained  below  show  that  the  interaction 
potential  between  gas  and  the  substrate  atoms  needs  further  careful  investigation, 
specifically  in  connection  with  rather  general  classes  of  substrate  solids  which 
have  ionic,  metallic  or  other  bonding  characteristics. 

Using  the  canonical  partition  function  one  obtains,  as  usual,  an  expression 
for  the  chemical  potential  ui  of  Che  i-th  component 


„  -  -kT  /  31n  2  1 

1  \  3Hi  j  T.V.Nj 

(1) 

where  Z  -  e'F/M 

(2) 

-  £  e-*v'kT 

(classical) 

(3) 

V 

Z  being  the  partition  function,  F  the  Helmholtz  free  energy,  T  the  temperature, 

V  the  volume,  Nj  the  number  of  molecules  of  the  j-th  species  (component) ;  E^ 
refers  to  the  energy  eigenvalue  of  the  v-th  state  associated  with  the  Hamiltonian 
H  of  the  appropriate  Schrodlnger  equation  of  the  system  in  the  quantum  mechanical 
partition  function 


4  being  the  matber  of  degress  of  freedom  of  the  system. 


(4) 
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For  Che  present  problem,  Che  potential  energy  U  In  the  Hamiltonian 
3  Pi2 

H  -  +  a<V***»V 


(5) 


is  the  total  effective  interaction  potential  between  the  N  gas  molecules  and  the 
solid  substrate  taken  as  a  whole. 

The  passage  from  quantum  statistical  mechanics  to  classical  statistical  mechanics 
referred  to  above  in  terms  of  the  respective  partition  functions  leads  to  a 
corresponding  transition  in  the  Helmholtz  free  energy 


F 


quantum  mechanical  classical 

thereby  forcing  the  following  factorization 


(6) 


-F  ,  /kT 

class. 

-Fideal/kT  e_Fconf  ig/lcT 

(7) 

e  * 

e  • 

with  Z  , 

class. 

Z ideal  ^conflg. 

(8) 

where  Zconfig. 

* 

V 

(9) 

and  Q  3  f... 

/  -U/kT  .  . 

’I  e  dq2*.  dqH  . 

(10) 

Here  Q  is  the  usual  "configurational  integral",  the  kinetic  energy  part  of  the 
Hamiltonian  having  gone  into  Zjjggj* 

We  know  (ref.  13)  that  the  chemical  potential  of  the  non- ideal  systems  breaks 
up  as 


“  *  “ideal  +  “ 


(11) 


where  y  is  the  part  arising  from  interactions  that  do  not  appear  in  the  ideal 
system;  thus 


9  In  Z 


ideal 


3  In  Z 


kT 


3N, 


/T.V.N 


3N, 


conflg. 


j 


/T.V.Nj 


(12) 


Since  we  are  interested  only  in  the  deviation  of  the  gas  adatoms  from  ideal  be¬ 
havior  and  beyond  that  in  the  variation  of  this  deviation  with  applied  stress — 
which  would  arise  from  a  corresponding  change  in  the  configurational  integral— 
we  need  only 


r 
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Zconflg. j 
Mi  /T.V.Nj 

la  (Q/vN3 

3Ni  jt,v,n.  ’ 


we  note  thereby  the  origin  of  the  factor  V**  that  appears  in  this  expression. 

The  validity  of  this  passage  to  the  classical  Unit  rather  than  the  use  of  the 
"Slater  sum"  (ref.  7)  is  rather  doubtful  when  the  solid  substrate  is  compressed 
or  more  generally  when  stresses  are  applied  to  it.  This  is  so  since  the  band 
structure  of  the  solid  can  change  considerably  involving  changes  in  che  ground 
state  and  furthermore  excitations  to  higher  levels.  Higher  states  thus  may 
contribute  to  changes  in  the  polarizabilities  of  the  solid  which,  as  we  shall 
see,  can  alter  the  interaction  potential.  However,  in  harmony  with  our  goal 
which  is  to  show  that  a  chemostress  effect  may  exist  based  on  general  grounds, 
we  shall  use  this  classical  limit. 

The  interaction  potential  U  is  usually  regarded  as  being  made  up  of  two-body 
interactions  between  any  adatom  and  each  of  the  atoms  of  the  solid  substrate. 
Furthermore  it  Is  assumed  that  the  total  potential  is  additlvely  (ref.  14)  made 
up  from  these.  While  this  assumption  may  yield  a  fair  approximation  to  the  actual 
interaction  for  extreme  ionic  solids,  it  certainly  falls  for  metal  substrates. 

For  these  one  needs  to  use  image  potentials  which  are  further  corrected  (ref.  15) 
for  effects  which  become  rather  significant  at  close  distances  between  the  adatom 
and  the  metal  substrate.  In  any  case,  there  will  be  an  effective  interatomic 
potential  u^  between  the  1-th  adatom  and  the  solid  stubstrate  taken  as  a  whole. 
Denoting  by  u^  the  effective  interaction  between  any  two  gas  molecules  and 
neglecting  higher  order  terms  we  have 

N  N  N-l 

U  «  £  u,  +  2  2  u,.  .  (15) 

i-1  1  j>k  k-1  J 

As  usual,  (ref.  16)  one  defines  the  functions: 

-u  /kT 

ft  =  e  1  -1  (16) 

-u. ./kT 


in  terms  of  which  the  configurational  part  of  the  partition  function  becomes 


(18) 


(19) 

(20) 


Here  we  pause  to  examine  the  meaning  of  the  volume  V  that  appears  above. 

As  discussed  by  Steele  (ref.  7),  the  volume  V  can  be  divided  into  two  parts;  a 
part  Va  which  is  a  small  volume  near  the  surface  where  the  gas  density  is  appre¬ 
ciably  different  from  the  density  in  the  remainder  (ref.  7)  V-Va  of  the  volume  V; 
then  the  main  contribution  of  Eq.  18  comes  from  Vft.  Denoting  the  totality  of  the 
surface  of  the  solid  by  A, 


X 


V  -  A.h  (21) 

a 

where  h  is  a  derived  distance  beyond  which  the  interatomic  potential  at  the 

temperature  T  of  the  system  has  no  effect.  Another  way  of  stating  this  1b  that  j  > 

m0 

beyond  a  perpendicular  distance  h  away  from  surface,  the  attractive  part  of  the 
potential  Is  significantly  less  than  kT  below  the  sero  of  the  potential  that  was 
taken  to  be  atinf inity»_gpd  therefore  these  molecules  are  essentially  free  from 
the  interaction. 

Another  parameter  that  is  significant  is  the  distance  zq  of  closest  approach 


between  adatoms  and  the  solid  atoms.  This  is  usually 
where  the  repulsive  part  of  the  potential  intersects  i 
dlcular  to  the  surface. 

Neglecting,  for  the  present  analysis,  the  and 
we  obtain 


Using  the  experimentally  and  theoretically  determined 


Eq .  22  becomes 


taken  to  be  approximately 
:he  z  axis  that  is  perpen- 

higher  order  contributions. 
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and  we  see  in  Fig.  1  that  y  is  sensitively  dependent  on  the  magnitude  of  the 

It 

minimum  of  the  potential  denoted  by  e  which  is  the  heat  of  adsorption.  We  see 

* 

from  Fig.  2  that  the  dependence  of  n  on  e  is  insensitive  to  the  powers  of  inverse- 

* 

power-law  interaction  potentials  and  that  log  n  is  linearly  dependent  on  e  .  One 

may  again  question  the  validity  of  such  general  conclusions  but  we  accept  it  since 

we  are  seeking  an  "average  thermodynamic"  effect. 

From  Fig.  10-1  of  Reference  7  we  observe  that  h  is  a  small  multiple  of  zq  and 

therefore  y ,  or  expf-y/kTl  which  is  a  more  significant  quantity,  becomes  sensitive 

* 

to  variations  of  the  depth  e  of  the  potential. 


THE  CHEMOSTRESS  EFFECT 

In  this  section  we  show  that  stresses  applied  to  the  solid  may  change  the 
chemical  potential. 

A  variety  of  interatomic  potentials  proposed  for  physical  adsorption  have  been 
reviewed  by  Margenau  (ref.  14).  For  extreme  ionic  solids,  additivity  of  the  two- 
body  interaction  between  adatom  and  substrate  atoms  can  be  assumed  to  be  approxi¬ 
mately  valid,  giving  the  result 


a  s  ^ 


(25) 


+  E 


where  D  is  the  perpendicular  distance  between  the  adatom  and  the  solid  surface, 
as  and  aa  are  the  polarizabilities  of  the  solid  and  of  the  adatom  respectively; 

E  is  the  mean  excitation  energies  of  the  adatom  and  Eg  is  that  of  the  solid. 
Applied  stresses  can  change  the  polarizabilities  of  the  solid  (ref .  17) .  In  the 
range  of  3-10  kilobars,  the  empirical  Jarman  rule  (ref.  IS) 


i  (t)T'2-°  <2' 

seems  to  hold  for  alkali  halides  in  general,  where  a  is  the  polarizability  and 
V  the  volume  of  the  solid.  This  translates  to 


~  %  2.0  3  (AP) 


(27) 


where  dP  is  the  pressure  change  and  g  is  the  isothermal  compressibility 


For  example,  for  NaCl,  2X  change  in  volume  for  vlO  kilobar  increase  in  pressure 
appears  to  be  (ref.  19)  a  reasonable  order  of  magnitude  estimate  based  on  experi¬ 
mental  data.  Therefore  a  change  of 
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%  0.04 


may  result  In  the  polarizability.  Referring  to  Eq.  25  above,  if  one  neglects 

variations  of  Eg  with  pressure  for  the  moment  (to  be  discussed  below),  a  fractional 

change  of  4%  in  the  polarizability  translates  into  a  similar  change  in  the  depth 

of  the  potential,  i.e.,  the  heat  of  adsorption  in  this  case.  The  experimental 

findings  of  Sams  et  al  (ref. 10)  then  clearly  indicate  that  that  ln(B._/Az  )  varies 
^  Ao  O 

linearly  with  e  /RT  giving  approximately  a  20%  change  in  B^/Az^,  if  we  assume 
that  these  experimental  findings  of  Sams  et  al.  can  be  extended  to  alkali  halides; 
this  is  a  rather  strong  assumption  that  is  partly  based  on  the  insensitivity, 
mentioned  above,  of  this  dependence  on  the  exponents  of  the  Lennard-Jones  type 
potential  as  well  as  on  the  universality  of  a  law  suggested  by  the  results  shown 
in  Fig.  1. 

Combining  these  results,  we  have 

_  Ajj.  ^  V**  in  (30) 

Ao  ^  1  +  (z  /ti)n  3a 


where  a  stands  for  a  general  form  of  stress  whose  tensorial  aspects  are  here 
neglected  in  line  with  the  purpose  of  the  present  paper.  Depending  on  the  values 
chosen,  this  then  may  give  a  1-10%  effect. 

There  remains  the  question  of  the  variation  of  the  mean  excitation  energy  Eg 
of  the  solid  with  applied  stress.  We  do  not  know  at  this  stage  whether  this 
change  (which  can  exist  due  to  changes  in  the  band  structure  when  the  solid  is 
stressed)  would  increase  or  decrease  the  magnitude  of  the  chemostress  effect.  In 
any  case  we  know  that  these  two  variations  may  be  intertwined  and  at  high  pressure 
one  may  not  be  allowed  to  assume 


3U  3U  3a  .  3U  s 

3^  *  3i  3?  +  3T  '  (3l) 

s 

namely,  the  separability  of  the  two  contributions.  This  brings  us  back  to  the 
question  of  the  validity  of  the  pairwise  additivity  of  the  interatomic  potential. 
Considering  even  the  possibility  of  insulator-metal  type  transitions  under  high 
pressures,  we  want  to  go  from  solids  with  ionic  bonding  to  metals  and  see  whether 
in  this  latter  extreme  the  chemostress  effect  may  also  be  observable.  We  use 
Bardeen's  (ref.  15)  result  for  the  van  der  Waals  Interaction  between  a  neutral 
molecule  and  a  metallic  surface. 

0  .  .  £ll  -  (32 

1 2D  Ba+(Cez/2^) 

where  rg  is  the  radius  of  a  sphere  whose  volume  is  equal  to  the  volume  occupied 
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per  electron  In  the  metal.  A  change  In  volume  then  translates  into  a  change  in 

the  van  der  Haals  attractive  part  of  the  interaction  potential  and  therefore  in 

i  the  minimum  of  the  total  potential  which  includes  a  repulsive  part.  Although  C 

2 

itself  depends  on  rg,  one  may  assume  this  to  be  a  weak  dependence  and  that  Ce  / 2rg 

is  of  the  order  of  magnitude  unity.  Thus  a  few  percent  change  in  r  gives  a 
*  8 
similar  change  in  e  to  within  a  factor  of  two.  It  appears  then  that  a  chemostress 

effect  may  be  expected  in  most  solids  on  general  grounds. 

It  must  be  emphasized,  however,  that  in  a  real  solid  under  stress,  dislocations 
are  generated  whose  dynamics  (such  as  migration  to  the  surface)  should  alter  the 
effective  "average"  dispersion  forces.  It  would  require  rather  extensive  studies 
to  separate  out  the  contributions  to  changes  of  heat  of  adsorption  that  come 
Indirectly  from  such  imperfections.  However,  no  matter  what  the  detailed  mechanisms 
and  the  particular  pathways  may  be.  If  the  heat  of  adsoprtlon  is  changed  under 
stress,  then  there  will  be  a  corresponding  change  In  the  chemical  potential,  indi¬ 
cating  that  the  chemostress  coefficient  may  be  a  non-negligable  entity. 


«*/RT 


Fig.  1.  Fit  of  the  experimental  results  of  the  3-12  potential  model 
(see  reference  10). 
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Fig.  2.  Theoretical  dependence  of  n  on  e  /RT  for  various  power-law  interaction 
functions  (see  references  10  and  7). 
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